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Abstract 

We consider the general Wigner function for a particle confined to a finite interval and 
subject to Dirichlet boundary conditions. We derive the boundary corrections to the "star- 
genvalue" equation and to the time evolution equation. These corrections can be cast in the 
form of a boundary potential contributing to the total Hamiltonian which together with a 
subsidiary boundary condition is responsible for the discretization of the energy levels. We 
show that a completely analogous formulation (in terms of boundary potentials) is also possible 
in standard operator quantum mechanics and that the Wigner and the operator formulations 
are also in one-to-one correspondence in the confined case. In particular, we extend Baker's 
converse construction to bounded systems. Finally, we elaborate on the applications of the 
formalism to the subject of Wigner trajectories, namely in the context of collision processes 
and quantum systems displaying chaotic behavior in the classical limit. 



1 Introduction 

The Moyal- Wigner- Weyl quantization constitutes a phase space quantization method alternative 
to canonical or path integral quantizations, |p|-||13||. Instead of wave functions and operators one 
deals with quasi-distribution (Wigner) functions and ordinary c-functions in phase-space. The 
Weyl-symbol associated with a general operator A(x,p) is given by, |1|, [|, |12|| : 



A w (x,p) = — dU dr ] Tr{A(x,p)e i ^ +ir 'P}e-^ x - iriP . (1) 

In this framework the average of A(x, p; t) is evaluated according to the formula: < ip\ A(x, p; t) \ip >= 
J_™ dx dpA w (x,p;0)F w (x,p;t), where F w (x,p;t) is proportional to the Weyl-symbol of the 
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quantum density matrix, p, and is known as the Wigner function of the system, 
state, p = \ip >< the Wigner function can be shown to be of the form: 



F w (x,p;t) 



1 



+00 



dye- 2i ™/ n r(x-y-M(x + y;t), 



For a pure 



(2) 



where ip(x; t) is the solution of Schrodinger's equation: ih^(x; t) = Htp(x; t). The Weyl "transform" 
establishes a biunivocal correspondence between the quantum algebra A of observables with stan- 
dard operator product • and quantum commutator [, ], on the one hand, and the "classical" algebra 
A defined over the classical phase space T*M with a "star-product" * and a Moyal sine-bracket 
[, ] M , on the other hand. The two latter operations are given by, | 



A w (x,p) * B w (x,p) = A w (x,p)ef J B w (x,p) = A 



1 w 



x,p 




x,p 




(3) 



A w (x,p),B w (x,p) 



w 1 



M 



in v 



II" 



*B W -B w *A 



w 



-A w (x,p) sin I - J J B w (x,p), (4) 



where J = 



One realizes immediately that in the limit h — > 0, the star-product and 



d_d_ _ d_d_ 

dx dp dp dx 

the Moyal bracket become the ordinary product of c- numbers and the Poisson bracket, respectively. 
It is also easy to check that the dynamics of the Wigner function is governed by the Moyal bracket 
(4): 

9F w 



dt 



x,p;t) = H (x,p),F (x,p;t) 



M 



(5) 



where H w (x,p) is the Weyl symbol of the quantum Hamiltonian, H(x,p). 

If the time-independent wave function happens to be an energy eigenstate with eigenvalue E, 

Hip(x) 



i.e. JrLip{x) = Eip(x), then |15, [16, [17], [T8| the Wigner function F^(x,p) satisfies an equivalent 
*-genvalue equation with identical *-genvalue: 



H w {x, p) * F™ (x, p) = F™ (x, p) * H w (x, p) = EFg (x, p) 



(6) 



More generally, to every eigenstate \a > with eigenvalue a of some operator A, there is one and 
only one associated stargenfunction F^(x,p), which is a solution of the stargenvalue equation: 
A w (x,p) * F^(x,p) = aF^(x,p) [|, [l^, [18). This stargenfunction allows for the evaluation of the 
probability of measuring the eigenvalue a: 



V(A = a) = J dx J dp F w (x,p)F^(x } p). 



(7) 



It is important to emphasize that, in most cases, the advantage in the Wigner approach does 
not reside in solving specific problems. Rather, it should be regarded as a means to grasp certain 
conceptual aspects of quantum mechanics and its connection with classical mechanics. Nevertheless, 
in some situations it could be a better starting point for finding solutions to specific problems. The 
Wigner formulation is a useful tool to derive kinetic equations in particular regimes (dilute gas, 
weakly interacting particles). In nonequilibrium statistical mechanics this is a familiar approach, 
In collision processes the Wigner methods are also a useful tool 0, S]. Contrary to the standard 
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operator approach they display the important advantage of making possible the use of approximation 
methods similar to the ones used in the full classical treatment. Furthermore the related subject 
of Wigner trajectories M, [|, 21]| provides a pictorial and practical tool to study collision process 



and it is of similar importance in the context of quantum chaotic systems [|], [U| [22|, £3|, |24[ . 

Boundary value problems, on the other hand, are one of the cornerstones of quantum mechanics. 
These are much more realistic models and moreover they entail a discretization of observables' 
spectra, which is one of the main features of quantum mechanics. They appear in virtually all 
branches of physics, ranging from quantum mechanics to general relativity, open strings and D- 
branes, where the non-commutative *-product and the Moyal bracket also find several applications 
[|16| , |25| , p6| . In standard operator quantum mechanics some famous examples of boundary value 
problems include the Kondo problem |27fl , quantum Hall liquids with constriction |28| and the 
Callan-Rubakov model . Furthermore one of the simplest collision process is the collision against 
an impenetrable boundary and some standard examples of quantum chaotic systems are bounded 

!§• 

In spite of the importance of boundary value problems in general and of these models in par- 
ticular, Wigner quantum mechanics does not provide a self contained and consistent formulation 
of confined systems. The aim of this paper is thus to provide such formulation. More precisely 
our aim is to generalize the two key equations of Wigner quantum mechanics (eq.(5) and eq.(6)) to 
the case when the original wave function is confined to an interval and satisfies Dirichlet boundary 
conditions. 

Let us then consider a general confined system. In standard operator quantum mechanics, we 
first solve the Schrodinger equation in the bulk (with Hamiltonian H), subsequently impose the 
boundary conditions that pin down the physical solutions and finally impose the wave function to 
be zero beyond the bulk. Alternatively, we may be able to provide a boundary potential Vd that 
confines the wave-function and then solve the unconstrained system, with extended Hamiltonian 
Hp = H + Vd and the original conditions at the boundaries. 

To obtain the solution of the problem in terms of Wigner functions, we just have to insert the 
solution obtained in the Schrodinger formulation into the definition of the Wigner function. This 
was the method followed in ref. ||. But the question that remains is: how do we solve the problem 
within the Wigner- Weyl formulation? 

The first approach cannot be easily translated into the Wigner framework because the highly 
nonlocal character of the Wigner function produces nontrivial effects of the boundary on the bulk 
part of the Wigner function, so that the stargenvalue equation has to be altered, as we shall see 
(sections 2-4). 

The second procedure seems to be the natural starting point to derive, through the Weyl map, 
the Wigner formulation for the bounded system. However this is not an easy task, mainly because 
the formulation of confined systems in terms of boundary potentials is already problematic at the 
standard operator level. In section 6 we will consider this approach and a) consistently solve the 
problem at the standard operator level and b) use the Weyl map to recover the bounded stargenvalue 
equation derived in section 3. 

This paper is organized as follows. In section 2, we discuss some aspects of the boundary condi- 
tions satisfied by the Wigner function. In sections 3,4 and 5 we evaluate the boundary corrections 
to the *-genvalue equation (6) and the time evolution equation (5) directly from the general for- 
mula of the confined wave function. Furthermore, we a) extend Baker's converse construction to 
bounded systems and b) argue that the results obtained are compatible with the preservation of 
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the Wigner function's normalization under time evolution. In section 6 we develop the boundary 
potential approach for standard operator quantum mechanics and show that through the Weyl map 
this formulation yields the proper Wigner formulation of confined systems. In section 7 we discuss 
some practical applications of the formalism: the computation of Wigner trajectories for confined 
chaotic systems and for collision processes. Finally in section 8 we present our conclusions. 

2 Wigner function in a finite interval 

Let us assume that the domain of the wave function is the finite interval a < x < b with a < b, 
so that the wave function vanishes outside this interval. As we now argue, quantum confined 
systems cannot be treated, in the Wigner formalism, through the standard procedure of solving a 
differential equation and imposing the boundary conditions thereafter. In particular, the energy 
stargenfunctions for these systems do not satisfy the stargenvalue equation (6) with the Hamiltonian 
of the bulk. 

As an example, consider the free particle in the infinite well, with a = — b = —L/2. The 
stargenvalue equation would be of the form: 

^*F w (x,p)=EF w (x,p). 
2m 

If we solve the corresponding Schrodinger equation with Dirichlet boundary conditions we get for 
the fundamental state: 

-2^dx~^ x) = m%) MX) = U COS [lV ■ 

If we now determine the associated Wigner function using (2) we get H] that F (x,p) is given by 
eq.(39,40). We can easily check that F w (x,p) does not satisfy the previous stargenvalue equation 
and thus the only possible conclusion is that this equation has to be altered. Similarly, we can check 
that the time-evolution of the system does not satisfy the Moyal equation with bulk Hamiltonian. 

Let us start by studying the general properties of the Wigner function for a confined system. 
Since the wave function vanishes outside the interval }a, b[, it satisfies the condition: 

ip*{x — y)ip(x + y) = unless a < x — y < b and a < x + y < b. 

This means that: 1) for a < x < Xq = s ^-, we have: 

tp*(x — y)ip(x + y) = unless a — x < y < x — a, 
similary, 2) for xq < x < b: 

if)*(x — y)ip(x + y) = unless x — b < y < b — x, 
and finally, 3) for x < a or x > b: 

ip*(x-y)tp(x + y) = 0, \/y. 
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We conclude that the Wigner function (2) for the bounded system can be written as: 



1 r+oo 

F w (x,p) = — / dye~ 2tpy / h r(x - y)ip(x + y) = G^F^fap) + e 2 (x)F™(x,p), 

Till J-oo 



(8) 



where: 

1, if a < x < x 

0i (x) = 9(x - a) - 9(x - x ) — < 

0, otherwise 

J 

1, if xo < x < b 

e 2 (x) = 6(x - x ) + 9(b - x) - 1 = < 

0, otherwise 

9{x) is the Heaviside's step function: 9{x) — 1 iff x > and 9{x) = otherwise, and finally: 

' F?(x,p) = ^j:: x a dye-^r(x-yMx + y), 

< 

k F™(x,p) = ± t-l dye- 2i ™/ n r(x - y)4>(x + y) 



(9) 



(10) 



Since the wave function or its derivatives are possibly discontinuous at x = a and x — b, and to 
avoid possible misinterpretations in our future calculations, the previous integrals are defined as 
improper: J^J" stands for lim c _ +(l + and likewise J^Zb stands for lim c _, fe - ^-c- Notice that this 
is fully compatible with our previous results identifying the domain where ip*(x — y)ip(x + y) is not 
identically zero. 

As a second remark let us point out that the boundary conditions on ip(x) are imposed on the 
bulk side of the boundary. For instance, the Dirichlet boundary conditions would be of the form: 
lim e ^ + if; (a + e) = lim e ^ + ip{b — e) = 0, which is not the same as requiring if) (a) = ip(b) = 0, 
since these are valid independently of the boundary conditions satisfied by if)(x). This is so because 
if)(x) is confined to the open interval ]a, b[. It is important to realize that we could equally choose 
to defined "confined" as "confined to the close interval [a, b}" and in fact for Dirichlet boundary 
conditions the two prescriptions yield exactly the same wave function. However, the confinement to 
an open interval makes some of the future steps more natural (though occasionally more involved) 
and this is why we choose to work with this prescription. 

Finally, and to make our future expressions more compact, we introduce the notation: /(c ± ) = 
lim e ^o+ f{ c =t e ) which will be used whenever there is no risk of misunderstanding. 

Let us then study the boundary conditions that are satisfied by the Wigner function of a confined 
system. From (8-10), one realizes immediately that: 



' F w (a+,p) = FW( a + lP ) = F^(b~,p) = F w (b',p) = 0, 
, F w (x^p) = F^(x ,p) = Ff (x ,p) = F w (x+,p), 



(11) 



where, accordingly to the previous notation, F w (a + ,p) stands for the limit lim e ^ + F (a+e,p) and 
likewise for the other expressions. Since F w (x,p) is defined as an improper integral it already re- 
quires the evaluation of a limit. This limit will always be calculate before the limit \im x ^ a + F w (x, p), 
i.e. first we evaluate the Wigner function for all x and only then do we compute whatever limits of 
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the Wigner function. Let us proceed: eq.(ll) means that F w (x,p) is continuous at x and obeys 
Dirichlet boundary conditions irrespective of the particular boundary conditions satisfied by the 
associated wave function ijj. The boundary conditions obeyed by the confined wave function act 
only on the derivatives of the Wigner function. For a < x < x$ a straightforward computation leads 
to: 

^(x,p) = ± {e-^- a V*(a + M2x - a+) + e~^ a - x ^*{2x - a + )^(a+)+ 



+ dye- 2 ^ h W*(x - y)il>(x + y) + r(x- yW(x + y)}} , 
where the prime denotes the derivative with respect to the argument. Consequently: 



\im x ^ a+ ^(x,p) = ± 



dx 



J_ 



e r 



(b-a 



)^*(a+)V'(6") + e~^ a - b ^*{b-)^ 



+ 



(12) 



(13) 



+0+ Ja-b 
2 



Similarly, from we get: 



liuu 



9x 



+ 



(14) 



+ ^ lim e _ >0+ rfye- 2 ^/ ft [^'*(x - y)^(x + y) + ^*(x - yW(x + y)\ . 

2 ~ l ~ 

We conclude that in general ^^-(xq >p) ^j^~( x o,p), exce Pt f° r Dirichlet boundary conditions 
(i.e. ip(a + ) = ip*(a + ) = tp(b~) = ip*(b~) = 0), in which case the Wigner function also obeys 
Neumann boundary conditions: ^-(a + ,p) = ^j^-(b~,p) = 0. We can carry on with this process 
and compute the second and third derivatives. The results concerning the boundary conditions 
can then be summarized as follows. All Wigner functions for bounded systems satisfy eq.(ll). 
These are consistency conditions and, given their general character, cannot give countenance to 
the discretization of the energy spectrum. The boundary conditions satisfied by the confined wave 
function yield subsidiary conditions on the Wigner function: 



Dirichlet: ip(a + ) = ip*(a~* 



ijj(b-) = ip*(b 
dF w 



a ,P) 



- 

dF w 



(b-, P ) = o. 



dx ' dx 

If a Wigner function satisfies these Neumann conditions, then it automatically verifies: 



(15) 



8F W 

dx 



x ,p) 



8F 



^z-(x ,p) 



Qx V^O 

d 2 F w 
dx 2 

Q2 F W 



(b-, P ) = 0, 
&p-( x o,P)- 



(16) 



Obviously, this pattern does not continue indefinitely. For instance, (a + ,£>) = -^%W{ a+ )\ 2 7^ 0. 
It is also important to emphasize that contrary to eq.(15), the additional set of eqs.(16), being just 
a consequence of eq.(15), do not constrain F w any further. 



Neumann: ^'(a+) = i/j'*(a+) = ^'{b') = ^'*(b~) = 



8F l 



dx 



(a + ,p)^0, 



dF l 



dx 



^(, -,p)^^(4,p) 

^(a+p) = SgTl(6-p)=0 



(17) 



fe 2 



(z , P) 



04, p). 



We also have: (a + ,p) ^ 0. 



The former results show that the confinement of the system and the boundary conditions satisfied 
by the wave function influence the Wigner function in a very non-trivial way. Furthermore, the 
conditions (11,15,16) will play an important part in the derivation of the new *-genvalue equation 
for bounded systems. 

However, one should realize that the previous results also display some less interesting features: 
let us focus on the case in which the wave function satisfies Dirichlet boundary conditions. If the 
system is also confined then the Wigner function will satisfy the set of conditions (11,15,16). The 
problem is that the converse result is not valid, i.e. the conditions (11,15,16) do not imply the 
confinement of the system and it is also not clear if, in general, they imply that the original wave 
function satisfy Dirichlet boundary conditions. Hence, (11,15,16) will not completely pin down the 
Wigner transform of the confined wave function satisfying Dirichlet boundary conditions and thus 
eqs. (11,15,16) fail to provide the proper translation (into the Wigner language) of the Dirichlet 
boundary conditions on the wave function. 

We now introduce a new set of boundary conditions for the Wigner function, which are fully 
equivalent to Dirichlet boundary conditions on the original wave function, both for the confined as 
well as for the unconfined case. We call these conditions integral Dirichlet boundary conditions and 
they can be derived as follows: if the wave function satisfies Dirichlet boundary conditions then 
lim e ^ 0+ V[a + e) = lim e ^ + \^{ a + e ) | 2 — and similarly lim e _ >0+ V(b — e) — lim e ^ + \i>(b — e) | 2 = 
where V(x) stands for the probability distribution. These equations can be written in the Wigner 
context quite straightforwardly: 



' \im^ Q+ J^dpF w (a + e,p) = 0, 
k hm^ S£dpF w (b-e,p)=0. 



(18) 



Notice that the use of the limits is valid (though unnecessary) in the unconfined case (where the 
wave function and its derivatives are continuous) but it is necessary if the system is confined. 

From the two sets of boundary conditions satisfied by the Wigner function, (eqs. (11,15,16) 
and eqs. (18)), the Dirichlet integral ones are the most appealing both from the physical and the 
mathematical points of view: firstly, because they represent a physically meaningful imposition, 
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that the probability of finding the particle at the boundary is zero. Secondly, because they are fully 
equivalent to the Dirichlet boundary conditions on the wave function. From eq.(7,18), we get: 

/+oo 
dpF w (a + e,p) = =>• lim \ip(a + e)| 2 = =>• lim + V(a + e) = 0, (19) 

and similarly for the 6-boundary. Finally, and to avoid a possible misinterpretation let us point out 
that eq.(ll) does not imply the satisfaction of the conditions (18) given the fact that, in general, 
it is not possible to interchange the order in which one evaluates the limit and the integral in (18). 
For instance, if a confined wave function satisfies Neumann conditions then the Wigner function 
will satisfy (11) but not (18). 



3 Boundary stargenvalue equation 

Let us henceforth assume that the wave function satisfies Dirichlet boundary conditions, so that 
eqs. (11,15,16,18) are valid. Furthermore, we assume that ip(x) is an energy-eigenstate for a Hamil- 
tonian of the form: H(x,p) = 2- + V(x). If we substitute F w (8) on the left-hand side of equation 
(6), we get: 

h w *f w = h w * {e^? + e 2 Ff } = { t- + v(x) \ * (e^r) + { £- + v(x) } * (e 2 Ff ) . (20) 



(21) 



(22) 



A straightforward calculation leads to: 

H w * ) = 0x (H w * F™) - [5(x -a)- 8(x - x )} F™+ 

+ £ [8'(x - xo) - 5'(x - a)] F™ + £ [5(x - x ) - 8(x - a)] , 
taking into account the boundary conditions at x — a, we get: 

H w * (0^) = X (H w * F?) + it5(x - x )F^(x ,p) + 

- *o)^(*o,p) + - x )FF(x,p). 

Similarly, we also have: 

H w * (0 2 Ff ) = 2 (H w * Ff ) - ^5(x - x )F™(x ,p)- 

- *o) 9J ^(x ,p) ~ - x )F?(x,p). 
Adding up all contributions and taking into account the fusing conditions (11) and (16), we obtain: 

H W *F W = 0! (H w * F?) + 2 (H w * Ff ) . (24) 

Let us then calculate the term H w * F^ for a < x < xq and the term H w * F^ for xq < x < b. 
The following theorem will do this. 



(23) 
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Theorem 

Let and be given by eq.(10) where ip{x) is an energy eigenstate satisfying Dirichlet boundary 
conditions at a and b. Let H w = j — h V(x). We then have: 



2m 
->W( 



H w *F^(x,p) = EF^(x,p) + ^5'(x-a)*F^(x + ,p) } for a < x < x 
H w *F^(x,p) = EF^(x,p)-^5'(x-b)*F^(x- J p) J for x < x < b 



(25) 



Proof 



Let us consider the product if * f\ . Following the method described in ]T5j] , we get: 



# w # F w = i 

J- 7r7i 



pa;— a j 

L-x dye 



ip*(x — y)ip(x + y) 



i 



2rmrh 



± J^dye-^V(x + y)^*{x - y)^(x + y) + 



p 2 - ihpd x -^d x 



(26) 



fi-2dye— ! Pil>*(x - y)ip(x + y). 
where d x acts on F^ and d x acts on H w . Taking into account that pe ^ = Y~§y e ^ an< ^ ^ na ^ 



ip(x) is an energy eigenstate, we get after a few integrations by parts: 

H w (x,p) * F^(x,p) = EF™(x,p) + B 1 {x,p). (27) 
The extra term Bi(x,p) is a boundary correction given by: 



Bi(x,p) 



2/p 



2nm 



(a 



i)*(2x - a~)i){a + ) + ip'*(2x - a")^(a + ) + tp*(2x - a~)ip'(a + ) I . 



For Dirichlet boundary conditions, we get: 



27rm 



^ (a -"V*(2x-a-)^(a + ). 



(28) 



(29) 



Let us now attempt to express B±(x,p) in terms of the Wigner function F^(x,p). Consider the 
following integral (where e > 0): 

A e = f+™ dke^-^kF™ (x + e,p - fife) = 



± dA; Sa-x-e dye ik ^ x - a+ y^ke- 2i Py^*{x + e - y)if>(x + e + y) = 
-- | ! x a Z:ll dy [p{x -a + y)] e^VO* + e - y)^(x + e + y) = 
-| S:i£: dy6(x - a + y) e -«W/» {-f + * - vM* + e + y)- 
+ e - y)^^ + e + y) + ^*(x + e - y)i>'(x + e + y)} = 



= f e-TM {_^*(2x - a + e)^(a + e) - ^'*{2x - a + e)^(a + e) + V*(2x - a + e)^'(a + e)} . 

(30) 
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In particular for a wave function satisfying a Dirichlet boundary condition, we conclude that: 
lim e ^ + Af (x,p; a) = ^e~^ ( - a ~ x ^ip*(2x — a~)ip'(a + ). It is straightforward to obtain from eq.(29): 
B®(x,p) = lim e ^ + ^^-Af (x,p; a). Let us now try to cast this expression in terms of ^-products: 

(31) 

= -2»E: = o h. (f y$ (n+1) (x - a)^f(x + c,p) = -2m [S'(x - a) * F^(x + e,p) . 

Taking the limit e — > + of the previous expression and substituting it in (29) and (27) yields 
the first equation of (25) directly. As a final remark we notice that 8'(x — a) * F^ v (x + ,p) = 
\im e ^ {5'(x—a)*F^ v (x+e,p)} is not identical to S'(x—a)*{\im t ^ F^ v (x+e,p)} = S'(x— a)*F^ (x,p). 
This last expression is not even well defined (cf.(30)) since the Dirac delta function is evaluated at 
a point where ip\x) is not continuous. Hence, we should keep in mind that in general, we cannot 
interchange the order in which the starproduct and the limit e — > are evaluated. 

Finally, we perform the analogous procedure for F 2 W and get the second equation of (25), which 
concludes the proof. 

Let us then return to equation (24). Using eq.(25) we get: 

H w (x,p) * F w (x,p) = 0i (x) {EF^(x,p) + ^5'(x - a) * F™(x + ,p)} + 

+& 2 (x) {EF 2 w (x,p) - ^5>(x - b) * F?(x-,p)} . 



(32) 



(33) 



Let us try to find the corresponding compact version for F w (x,p) (eq.(8)). We notice that: 
1) QiEF^ + Q 2 EF 2 W = E {QiF^ + G 2 F 2 W } = EF W . 
2) 

0i (x) {S'(x -a)* F™(x + ,p)} = e 1 (x) {-±re-^ a - x ^*(2x - a")^'(a+)} = 

= &i(x + ) {-^e-^ a - x ^*(2x - a-)ip'(a + )} = 8'{x - a) * {0 1 (x+)F 1 w \x+ ,p)} , 

where in the first step we used eqs. (29,31); in the second step the fact that Q\(x) = lim e ^o+ &i(x + 
e), Vrr ^ a, Xo and that ip*{2x — a~) vanishes for x = a and x = x ; and finally in the last step that 
Qi(x + ) is p-independent. 

3) Similarly, we have: Q 2 (x) {S'(x - b) * F 2 w (x~,p)} = S'(x - b) * [Q 2 {x~)F 2 w \x~ ,p)}. 

4) Finally, we also have: 

5'{x - b) * {e 1 (x + )F] v (x + ,p)} = 0i(x + ) {5'(x -b)* F^(x+,p)} = 
= £;Qi(x + ) /+~ dke lk ^kF^ (x+,p-\k) = -^ K Q 1 (x + )6(x + -x )e- 2 -^^x (34) 

x {-^*(2x - b~)ij(b + ) - ^'*(2x - b~)^{b+) + ^*{2x - b-)^{b+)) = 0, 

where in the third step we made the same calculation as in eq.(30) and in the last step we used the 
fact that ip(b + ) = ip'{b + ) = 0, or else that Qi(x + ) = 0, Wx > x and 9(x + — x ) = 0, Vrr < x . 
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5) Similarly: 



0. 



6'{x-a) * (e 2 (x-)F^(x-,p)] 
Using these results we can re- write eq.(32) as: 

H w *F W = EF W + - a) * {6i(x+)Ff (x+,p) + e 2 (ar)Ff (aT,p)} - 

-£;5'(x - b) * {Q 1 (x + )F^(x + ,p) + e 2 (x-)F^(x-,p)} , 
from where it immediately follows: 



(35) 



(36) 



p 

2m 



+ V(x) *F w (x,p)-—5'(x-a)*F w (x + ,p) + —5'(x-b)*F w (x-,p) = EF w (x,p). (37) 



2m 



2m 



Following the same procedure, one would also obtain: 



F w {x,p)*(2- + V{x) 
\ 2m 



Ti 2 ft 2 

—F w (x + , p) * S'(x -a) + —F w (x~ , p) * 5'(x -b) = EF w (x, p). (38) 
2m 2m 



We conclude that, if the eigenfunction satisfies Dirichlet boundary conditions, then the correspond- 
ing Wigner function does not satisfy the standard *-genvalue equation (6), but rather a modified 
equation given by (37,38). Furthermore, the source of discretization of the energy spectrum are 
not the Dirichlet conditions on the Wigner function but instead the subsidiary boundary conditions 
(15) or (18). 

As a simple example let us consider the free particle in the infinite potential well. We choose as 
Ansatz: 



F w (x,p) = -ir p cos [2(P-k\x\)]-sm 



2£ ( k _ 
h \2 



x\ 



+ 



sin 



s(p + **)( 



\x\ 



+ 



4Tr(p—hk) 



sin 



(p - hk) (j 



(39) 



— \x 



corresponding to the wave-function if)(x) = asin(kx + (3). This function is continuous at xq = 
and satisfies Dirichlet 's condition at x = ±-|, in accordance with (11). If we impose the constraints 
(15) and substitute this expression in the *-genvalue equation (37,38), using the integral form A 
(eq.(30)), we obtain: 



kr. 



nir 



Pn 



nn 



2/„2 



En 



2^2 



h 7T 



2m 2mL 2 

rb J_ f + OO j„ J7i 



n 



71=1,2, 



(40) 



Finally, if we impose the normalization f a dx f_™ dpF (x,p) = 1, we obtain a = J2/L. This result 
is in perfect agreement with ref. [H]. 



4 Baker's converse construction 

In the last section we proved that if the wave function is confined to the interval ]a, b[, satisfies the 
operator eigenvalue equation and Dirichlet boundary conditions at a + and b~ then the corresponding 
Wigner function satisfies a modified *-genvalue equation together with integral Dirichlet boundary 
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conditions at a + and b~. In this section we want to prove the converse result. In the unconfined 
case this result is known as Baker's converse construction || HJ. Hence, in this section we want to 
extend Baker's converse construction to the confined case thus proving the full equivalence of the 
Wigner and operator formulations of quantum mechanics when boundaries are present. 
Let us then consider some real function F (x,p) satisfying the following conditions: 

(i) it is a continuous function of x and p; 

(ii) it has an infinite number of derivatives with respect to p and is twice differentiable with respect 
to X] 

(iii) F w obeys integral Dirichlet boundary conditions at x = a + and x = b~ , eq.(18); 

(iv) it obeys the left- and right- boundary stargenvalue equations (37,38) when e — > + : 

ii + V ( x )) * F w (x,p) - ^5'{x - a) * F w (x + e,p) + ^5'(x - b) * F w (x - e,p) = 



= F w (x,p) * (£ + V(x)) - ^F w (x + e,p)* 5'(x - a) + £f w (x - e,p) * 5'{x -b) = EF w (x,p). 

(41) 

We will now prove that under these assumptions, there is a unique (up to a global phase factor) nor- 
malized complex and continuous function tp, related to F w according to (2), which obeys Dirichlet 
boundary conditions at x = a + and x = b~ and is a solution of Schrodinger's equation. 
It will prove useful to consider the Fourier transform of F w : 



T r+oo 

F w (x,p) = — dy e- 2ipy / h F(x,y). 
irk J-oo 



(42) 



Substituting this expression in the left- and right- stargenvalue equations (41), and following the 
same steps as in eq.(26), we obtain: 



in dy e" | [_£ (^) 2 + V{x ±y) _ E 



F(x,y) 



(43) 



-^S'(x ±y- a)F(x + e,y) + ^8>{x ±y- b)F(x -e,y)} = 0. 
This means that F(x, y) satisfies the equations: 



h 2 ( d x ± d„ 



2m 



+ V(x ± y) 



h 2 ~ h 2 

F(x,y)-—5'(x±y-a)F(x+e,y)+—5'(x±y-b)F(x-e,y) = EF(x,y). 

(44) 



From now on we will follow the steps of the standard Baker converse construction for unconfined 
systems ||. Introducing the function F such that F(x — y,x + y) = F(x,y) and performing the 
change of variables u = x — y and v = x + y, we get: 



2m dv 2 



1,1 >J + V(v)] F(u, v) - £S'{v - a)F(u + e ,v + e) + £S'(v - b)F(u - e ,v - e) = EF{u, v) 
-&& + V(u)] F(u, v) - %-8>{u - a)F(u + e,v + e) + £-5>(u - b)F(u - e, v - e) = EF{u, v) 



(45) 
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from where it follows that F(u,v) = £(u)ip(v) and so F(x,y) = £(x — y)i/)(x + y). Furthermore, 
since F w (x,p) is real, eq.(42) implies that £ = ip* and so F(x,y) = ip*(x — y)tp(x + y), where the 
complex function ip obeys the equation: 

fi 2 h 2 h 1 

i)'\x) + V(xU(x) 5'(x - a)^(x + ) + — 5'(x - bU(x~) = EMx), (46) 

2m 2m 2m 

where we performed the limit e — > + . Finally, if F w (x,p) satisfies integral Dirichlet conditions at 
x = a + and x = b~ then the corresponding wave function satisfies ip(a + ) = ip{b~) = 0, eq.(19). 

In summary, we have shown that if F w satisfies the conditions (i) to (iv) above then F w is the 
Wigner transform (2) of a complex function ip which 1) is a solution of the Schrodinger equation 
(46) and 2) obeys Dirichlet boundary condition at x = a + and x = b~ . In section 6 we will show 
how to solve this equation and prove that its unique (up to a global phase factor) solution is the 
standard solution of the corresponding eigenvalue problem in operator quantum mechanics. 



5 Boundary dynamics 

Let us now consider the time evolution: 

dF w 



dt 

and focus on the time derivative dF™ /dt: 



dt 



dt ' 



(47) 



dt 



^ dye 



-2ipy/h 



^r(x - y)ip(x + y)+ ip*(x - y)$(x + y) 



i 

ink 



i:: x a dye-^ h {£- [T(x - yWx + y)-r(x- y)f{x + y)] + 



(48) 



+ [V(x + y)- V(x - y)\ ^*(x - y)^{x + y)}, 

where we used the Schrodinger equation. Consider the first term. After a few integrations by parts 
and by keeping track of the boundary contributions and Dirichlet's condition, we obtain: 

^f:: x a dye-^r*(x-yMx + y) = ^{^ {x ~ a Vl« + )^(2^-«1- 

(49) 

-e-^ a -^*(2x - a-)4>>(a+)} + ^ dye~^r{x ~ vW'{x + y) - 



m dx 



The penultimate term exactly cancels the second term on the right-hand side of eq.(48). The first 
two terms in eq.(49) can be written as: 



h 



/-« dk 



Amir J-oo 



e ik(a-x) + e ifc(x-o)l k pW (x+ p _ h k 



ih 
'2m 



F^(x + ,p) * 5'(x - a) - 5'(x - a) * F^(x + ,p) 



(50) 



EL. 

2m 



5>(x-a),FW(x+,p) 



M 



The terms involving the potential V(x) are familiar in Wigner quantum mechanics. They yield the 
contribution 

dt 



V(x), F^ ¥ (x,p) . After assembling all the results, we get: 



(x,p;t) 



^ + V(x),FF(x,p;t) 



M 



2m 



5'(x-a),F^(x + ,p;t) 



M 



(51) 
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A similar calculation leads to: 



dFf_ 
dt 



x,p;t) 



^- + V(x),F^(x,p;t) 



+ 



M 



2m 



5'(x-b),F^(x-,p;t) 



M 



From these results we get for F w (x,p;t) = Q 1 (x)F{ v (x,p;t) + Q 2 (x)Ff(x,p;t) : 



dF 



w 



dt 



-(x,p;t) 



^- + V(x),F w (x,p;t) 
2m 



M 



+ 



2m 
2m 



5'(x-a),F w (x + ,p;t) 
5'{x-b),F w {x-,p;t) 



M 



M ' 



where we used the set of relations given by eqs. (33-36). 

A straightforward consequence is that probability is conserved (as expected): 



-j a dxj ^ dpF w (x,p;t) = -J ^ dx I dp Q 1 (x)F l w (x,p;t) + Q 2 (x)F 2 w (x,p;t) 



(52) 



(53) 



= 0. 



(54) 



6 Weyl transform of the standard quantum description 

The standard approach to derive the Wigner formulation of the eigenvalue problem for confined 
systems would amount to applying the Weyl transform, W eq.(l), to both the eigenvalue equation 
and to the boundary conditions. However, we saw in section 2, that this procedure does not yield 
the correct results. As we shall see, this is linked to the way operator quantum mechanics handles 
confined dynamical systems. The standard procedure is to 1) solve the eigenvalue equation for the 
unconfined wave function <f>(x) and 2) impose the wave function to be zero outside the bulk and 
produce the confined state ip(x) = <fi(x)9(x — a)9(b — x). The problem in the Wigner formulation 
is that W{\ip >< ip\) 7^ W(\<f> >< (j>\)0(x — a)6{b — x) and thus we cannot solve the unconfined 
stargenvalue equation and manipulate its solution to obtain the confined Wigner function. 

The reason for this apparent contradiction between the standard operator and Wigner formula- 
tions of quantum mechanics resides in the fact that ip is not a globally valid solution of the operator 
eigenvalue equation. Globally, ip satisfies a different eigenvalue equation which includes boundary 
potentials. Given the intrinsically non-local character of the Wigner function this fact is crucial to 
derive the correct Wigner description of the bounded system. 

In this section we will a) derive a globally valid eigenvalue equation for confined systems satis- 
fying Dirichlet boundary conditions, b) show that the only solution of this equation is the confined 
wave function and c) shortly discuss the numerical implementation of the new eigenvalue equation. 
It will then be trivial to realize that the Weyl transform of the new eigenvalue equation is the 
bounded stargenvalue equation (37,38), thus recovering the usual correspondence between standard 
operator and Wigner quantum mechanics. 



6.1 Bounded eigenvalue equation 

Let then be the unconfined solution of the eigenvalue equation satisfying Dirichlet boundary 
conditions at x = a and x = b: 

H(j> = E(j> , H=^- + V(x) and 0(a) = 0(6) = 0. (55) 
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9 e (x) 



Notice that <f>(x) is the solution of eq.(55) everywhere, i.e. for all x E 1Z. To obtain the confined 
wave function (for instance for a particle in a box) the standard procedure is to cut off (p outside 
the box and produce the new wave function: ip(x) = <f)(x)8(x — a)6(b — x). Notice that ip satisfies 
eq.(55) inside the box, it also satisfies it outside the box, but it does not satisfy it at the boundaries. 
Hence, ip is not a global solution of (55). Due to non-local effects the Weyl transform of ip does not 
satisfy the stargenvalue equation (6), even in the bulk. 

It seems quite natural to expect that if we are able to provide the global eigenvalue equation for 
ip then the correct stargenvalue equation would just be its Weyl transform. The aim of this section 
is to derive this global equation. The first step is to introduce a twice different iable version of the 
Heaviside step function and of the confined wave function ip: 

G C 2 (1Z) and is non decreasing 

= 0, x < — e , ipe( x ) — 4>(x)9 € (x — a)9 e (b — x). (56) 

= 1, x> e , 0<e<<l 

The definition of 9 e (x) is fully compatible with the distribution 6{x) we used in the previous sections: 
lim e ^ + 6e(x) = if x < 0, lim e ^ + e (x) = 1 if x > and lim e ^ + ^e(O) remains unspecified (it might 
be any number between and 1 or it may not be defined at all). The results of this section are not 
dependent of the particular function 8 e (x) but only of the general properties given in (56). From 
(56) the smooth versions of 5(x) and S'(x) follow immediately: S t (x) = Q' e {x) and S' t (x) = 0"(x). It 
is also clear that \im e _> +(5 t (x), S' e (x), ip € (x)) = (5(x),5'(x),ip(x)). Let us then apply H to ip e (x): 

HAM - {-££ + V{x) }m*) 

= {-^"(x) + V(x)<f>(x)} 0e(x - a)6 € (b - x) 

h 2 

- {2(P'(x)S £ (x -a)- 2(P'{x)6 e (b -x) + <p(x)5' e (x - a) + <f>(x)6' e (b - x)} 

= Erp e (x) + ^{<f>(x)6' e (x-a)-<j>(x)6' e (x-b)}, (57) 

where in the last step we use the fact that <p(x) satisfies the eigenvalue equation (55) and that the 
Dirichlet boundary conditions on (p{x) imply for sufficiently small e: <p'(x)5 e (x — a) = —<p(x)5' e (x — a) 
and likewise <p'(x)5 e (b — x) = —<p(x)5' t (x — b). The former equation is not yet a closed eigenvalue 
equation for ip € (x) since it also evolves the wave function <p{x). 

In what follows and to simplify the discussion we will consider the case of just one boundary 
placed at a = 0. The derivation of the two-boundary eigenvalue equation follows exactly the same 
steps and will be written explicitly at the end. In the one boundary case eq.(57) reduces to: 

HrP e (x) = Eip e (x) + ^-5' e (x)<f>(x), (58) 

where <p(x) satisfies eq.(55) with Dirichlet boundary conditions at x — 0. The former equation 
follows from applying H to the state ip e (x) = <f>(x)0 e (x). 
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To proceed we notice that for x > e the two equations (55,58) - for cf)(x) and for i/) e (x), respectively 
- are identical. If we supply identical boundary conditions, for instance, at x = 2e we will get 
identical solutions for x > e, i.e. if 4>(2e) = ij) e (2e) and (f>'(2e) = ip' € (2e) then (f>(x) = ip e {x) for all 
x > e. Furthermore, and since <f>(x) is analytical, we also have: 

00 (—2e) n 8 n 6 

0(x) = 0(x + 2e) + £ ^^(* + 26), (59) 

and for x > — e we have 0(x + 2e) = if) e (x + 2e) and equally §^(2; + 2e) = ^tt{x + 2e), Vn G A/". 
Hence: 

0(x) = ^(x + 2e) + £ ( ~ 2 ' r ^ (x + 26) , Vx>- e . (60) 
^ n\ ox n 

The trick here is that while <fi{x) is analytical ip e {x) is not (cf.(56)) and thus the Taylor expansion 
(60) yields (j>(x) instead of if) € (x). Substituting the former expansion in eq.(58) we get: 

ft 2 ( 00 (— 2e) n d n ib 1 

H^(x) = EA(x) + ^5' e {x) !^ e (x + 2e) + Y, 1 J ^ (x + 2e)| . (61) 

Notice that the previous equation is valid for all x G 1Z because S' e (x) = if x < — e. Finally by 
taking the limit e — > + we obtain: 

h 2 

Hij}(x) = EMx) + — 5\x)ij(x + ), (62) 
2m 

where we used the notation x + to make it explicit that in the product of 5' by ip the two factors 
are evaluated at different points. This is crucial because equation (61) is valid for all x G 1Z only in 
this case, i.e. in the limit e — > + but not for e = 0. Equation (62) is the new eigenvalue equation 
for ip(x). Following exactly the same procedure we can easily generalize it to the case of a double 
boundary placed at x = a and x = b, (a < b). We get: 

h 2 

HMx) = EMx) + — \5'(x - a)^(x + ) - 5'(x - bU(x-)} , (63) 
2m L > 

which together with the boundary conditions ip{a + ) = if)(b~) = provide the operator formulation 
of the confined eigenvalue problem for Dirichlet boundary conditions. To finish let us point out 
that: 1) eq.(63) was obtained in section 4 from the stargenvalue equations (37,38) and 2) conversely, 
the Weyl transform of eq.(63) immediately yields the stargenvalue equations (37,38). 



6.2 Solving the new eigenvalue equation 

Our next step is to solve the modified eigenvalue equation (63) and show that the confined wave 
function ip(x) is its only solution. To make it simple let us consider again the one boundary example 
eq.(62). 

For x > 0, 5(x) = and eq.(62) reduces to the unconfined equation (55). The boundary 
conditions ^(0 + ) = and ip'(0 + ) = 0'(O + ) impose ip(x) = <f>(x) for all x > 0. 
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For —a < x < a and a « 1 we have: 



J —a J — a Z l Jl J — a 

-^W {x) {-<?)} + 0{x + a) = 

EO{x + a) + ( ^(x')^' + + )l * - T ^(x')^ V + + )cfe' 



2m 

=^ ^'(x) - V'(-<t) =^'(0 + )^(x) + C(a; + ( T), (64) 

where in the first step we introduced the notation 0(y) to designate an arbitrary continuous function 
such that linij^o 0{y) = 0, and in the third step we used the fact that ip(0 + ) = 0. Taking the limits 
a — > + and x — > + we get: 

^'(0+) - V'(0") = V'(0 + ) =^ V'(0~) = 0. (65) 
The second integration of the eigenvalue equation yields: 

i/>'(x)dx- r ip'{-a)dx = f ip'{0 + )e{x)dx + f 0(x + a)dx 



=t> ip{(j) — ip{ — <t) = a{24>'(-a) +V'(0 + )} + O( ( r 2 ), (66) 

and in the limit a — > + we have -^(O" 1 ") — ^(O - ) — 0> which together with the original Dirichlet 
boundary conditions at x = + imply ^(0 _ ) = 0. 

We conclude that for x < 0, VK^) satisfies eq.(55) with boundary conditions ip{0~) = ip'(0~) = 
and thus ip(x) = 0, \/x < 0. Hence, we are left with the original bounded wave function: ip(x) = 
4>(x)9(x). The generalization for a system with two boundaries is straightforward. Following exactly 
the same procedure we solve eq.(63) and get ip(x) = <f>(x)9(x — a)6{b — x), as we should. This result 
concludes the proof of Backer's converse construction initiated in section 4. 



6.3 The free particle bounded eigenvalue equation: numerical solution 
of a simple example 

As a simple example let us consider the case of a free particle confined to the interval ] — 1, 1[ and 
subject to Dirichlet boundary conditions. 

The unconfined Hamiltonian in the position representation is given by H = -££r and the 
unconfined fundamental state (solution of eq.(55)) is: 4>i{x) = cos(|a; ) with E — Ei — h 2 n 2 /8m. 
Fig. 1 displays the unconfined wave function for h = m = 1. 

We now consider a first possible approximation to the confined eigenvalue equation (63). We 
make e = 0.25 and use a Gaussian approximation to the Dirac delta function: 5 e (x) = ^5727 exp | — (p-) 2 

where e' substitutes e for the following reason: the function S € (x) is not identically zero for x > e 
nor for x < — e and thus there is no function 9 e (x) satisfying the conditions (56) and such that 
S e (x) = 9' e (x). However, if we define 5 e (x) using a sufficiently small spread e' (when compared with 
e) then, for our numerical purposes 5 e (x) will provide a good enough approximation to 9' e (x). In 
all future applications we will make: e' = 0.8e. Furthermore, and to simplify the notation we will 
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refer to S e (x) as just 6 e (x). Finally, taking into account the considerable large value of e we have to 
consider the contribution of higher order terms in the expansion (60), and thus we will use the two 
boundaries generalization of eq.(61) as our starting point and truncate the series at the third term. 
Hence, the equation we want to solve is: 

-fifoMz) = E^{x) + + 1) {^ e (x + 2e) - 2ef £ (x + 2e) + 2e 2 <(x + 2e)} 

(67) 

-7^S' e {x - 1) {ij e (x - 2e) + 2ei[> e (x - 2e) + 2e 2 ?//'(x - 2e)} . 

Fig. 2 displays the numerical solution of the former equation for the boundary conditions tp € (a+2e = 
-0.5) = ip e (b - 2e = 0.5) = 0(0.5) = 0.707 and for E = E 1 = tt 2 /8, h = m = 1. Notice that the 
wave function is not completely confined which follows from the fact that we are not using the true 
Dirac delta but a (poor) Gaussian approximation of it. If we decrease the value of e the confinement 
will increase greatly. 

Let us then consider a second approximation to the confined eigenvalue equation by making 
e = 0.0025. Given the considerable small value of e we can now truncate the series in the two 
boundaries generalization of eq.(61) to zero order and numerically solve the finite version of eq.(63): 

' ' i/; e (x) = Eil> e (x) + — {6 , e (x+l)il> e (x + 2e)-6' e (x-l)il> e (x-2e)}. (68) 



2m dx 2 2m 

Fig. 3 displays the numerical solution of the former equation for the boundary conditions tp e (a+2e = 
-0.995) = ip e (b -2e = 0.995) = 0(0.995) = 0.008 and % = m = 1, E = E x = vr 2 /8. 

We will see in the next section that the part of the wave function that escapes from the box 
corresponds to quantum trajectories with energies much higher than the fundamental state. By 
making e — > + we completely confine the wave function. 



7 Wigner trajectories in confined systems 



In this section, we briefly review the topic of Wigner trajectories 0, [5|, [2(], |2T|], where the previous 
work may find some interesting applications. There are two well-known instances, where this concept 
plays an important role: collision processes M M and the study of quantum systems that exhibit 



chaotic behavior in the classical limit 0, PL3L 22, 23, 24 



The Schrodinger and Wigner formulations of a collision process are obviously equivalent. This 
notwithstanding, there are important practical differences between the two approaches. In most 
cases, where an exact solution to the problem is unknown, one has to devise some approximative 
scheme. If certain conditions hold (concerning the interaction potential), the Moyal equation (5) 
can be truncated to a good degree of accuracy. Since it does not contain any operator, it is often 
easier do develop approximations on it than on the Schrodinger equation. Having computed the 
Wigner function that solves this truncated version of the Moyal equation, one can subsequently 
evaluate quantum corrections to the classical phase space trajectories. In addition to the benefit 
of this pictorial and intuitive description in terms of trajectories there is another advantage in this 
approach: collision cross sections can be computed by applying the same numerical techniques used 
in classical collision processes. 

Another situation where the Wigner trajectories play an important role is that of nonlinear 
systems exhibiting classical chaotic behavior. In these systems the classical phase space trajectories 
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are of random nature while the Wigner trajectories are smooth paths thus showing that the chaotic 
behavior disappears upon quantization. A topic of current interest is that of understanding the 
transition from classical to quantum dynamics. Since the Wigner trajectories might be computed 
as an h order by order correction to the classical trajectories they provide a powerful, pictorial tool 
to study this transition. 

In this section we will a) shortly review the subject of Wigner trajectories for unconfined systems, 
b) discuss the main difficulties in extending the formalism to bounded systems and show that our 
previous results provide a new approach in this context and c) illustrate these results using the 
simple example of a free particle confined to an infinite square well. The reader should be advised 
that our aim is not to provide a fully developed and precise formulation of these topics but only to 
outline (and hopefully motivate) some possible lines of research. 

7.1 Wigner trajectories 

2 

Let us then consider the Moyal equation of motion (5) associated to H — ^ + V(x): 

dF w p dF w dV dF w (h/2i) 2 d 3 Vd 3 F w 

dt m dx dx dp 3! dx 3 dp 3 

Alternatively, we can write this equation as: 

dF w t ^ p dF w t /■+«> 
~dt 

where 



v Or r+oo 

( x ,p,t) = -^—-+ dp' J(x,p')F w (x,p + p',t), (70) 

771 OX J— oo 



J(x, p) = —2 / dy [V(x + y)- V(x - y)\ e~ 2ipy l h . (71) 

In the absence of a boundary and if the potential V is at most quadratic in x, then the Moyal 
equation (69) reduces to the Liouville equation: 

dF w p dF w dV dF w 

dt m dx dx dp 

The previous equation means that the solution F w (x,p,t) evolves along the classical trajectories, 

f F w (x,p,t) = F w (x(-t),p(-t),0), 



(73) 



This motivates the concept of Wigner trajectories. If we define an effective quantum potential 
V e ff(x,p,t) according to: 



dV e ff dF 



w 



dx dp 

then we can express the Moyal equation in the form: 



/+oo 
dp' F w (x,p + p',t)J(x,p'), (74) 



dF w pdF w | dV eff dF w 
dt m dx dx dp ' 
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which can be interpreted in analogy with (72), by stating that the Wigner function evolves along 
the Wigner trajectories given by: 



i= P (76) 

m dx 

If an exact solution for F w is known, then we can, in principle, calculate the effective potential 
from (74) and subsequently solve (76) to obtain the Wigner trajectories. However, in general, such 
is not the case. If the potential V(x) does not deviate appreciably from a quadratic potential then 
we can address the problem by iteratively evaluating order by order the quantum corrections to the 
classical trajectories j|, |2D|, |2"I| . The zero order approximation (i.e. the classical solution), satisfies 
the Liouville equation: 

d^ = _ L dF^ + d\^dFl_ 

dt m dx dx dp 

The order zero trajectories correspond to the classical Hamilton equations. The order one equation 
is given by: 

dF]" _ p dF^ | dV dF^ | (h/2if d 3 V d 3 F^ 
dt m dx dx dp 3! dx 3 dp 3 

If we solve this equation for F± and substitute in (75,76), we obtain the first order trajectories: 



x = — 

m ' 



(79) 

• 9V | 1 ( h\ 2 d 3 Vd 3 F™ ,dFf_ 

V dx " t " 3! V 2i ) dx 3 dp 3 I dp " 

The second order approximation is: 

dF™ _ p dFf + dV dF™ + (h/2if d 3 V d 3 F™ + (h/2i) 4 d 5 V d 5 F™ ^ 
dt m dx dx dp 3! dx 3 dp 3 5! dx 5 dp 5 ' 

and yields the second order correction to the classical trajectories: 



x = ^, 

m. * 



•_ dV | 1 / h \ 2 d s V 9 3 F^ i dFf i / a \ 4 py&F™ i dF); 
V dx "r 3! V 2i ) dx 3 dp 3 ' dp 5! \ 2i J dx 5 dp 5 ' dp 



51) 



This procedure can be carried on to any order in h. Its main advantage is that it transforms 
the problem of solving an infinite order partial differential equation (the Moyal equation) into 
the problem of solving a sequence of first order partial differential equations. It also displays 
the advantage of casting the transition from classical to quantum dynamics as an order-by-order 
correction (in H) to classical mechanics. There are, however, some weak points in the procedure, 
the most significant being perhaps that it is not always clear what is the degree of precision of an 
approximation to a given order nor if the entire scheme will converge. 
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7.2 Wigner trajectories in bounded systems 

If a system has a boundary, then the former approximative scheme will obviously break down, 
because the boundary contribution was thus far unknown. As we have seen in the previous sections 
the boundary interaction has a non-local nature and will thus affect the Wigner trajectories well 
inside the bulk part of the system. It therefore must equally contribute to the quantum corrections 
of the trajectories. 

In this case one was left with a single possible approach: to solve the Schrodinger equation 
for the confined wave function and compute the Wigner function from it. In some cases where an 
exact solution of the Schrodinger equation is known this method can be taken to completion. Take 
for instance the simple example of a free particle confined to the interval ] — 1, 1[ and subject to 
Dirichlet boundary conditions ||. The fundamental state is <pi(x) = cos(|x) from which we can 
compute the Wigner function eqs. (2,39,40). Since the particle is in an energy eigenstate we have: 

r)F w 

— = [H,F W ] M = =► F w (q(-t),p(-t),0) = F w (q,p,t) = F w (q,p,0), (82) 

and thus the Wigner function evolves along the paths where it displays a constant value, i.e. 
the Wigner trajectories (q(— t),p(— £)) are the equi- Wigner curves. Fig. 4 displays the Wigner 
trajectories for this simple system. 

For more elaborate examples we are forced to use eq.(74) or eq.(75) to calculate the effective 
potential from the Wigner function and then eq.(76) to calculate the Wigner trajectories. For 



instance in 0, ^4J a system in a symmetric box with Dirichlet boundary conditions and an external 
driving force was considered. The force, which is of the form F(t) = Focos(ut), induces classically 
the appearance of resonance. Since the system has a boundary, quantum mechanically, the authors 
were compelled to work in the context of the Schrodinger equation but since an exact, analytic 
solution is not known the Wigner function and subsequently the Wigner trajectories could not be 
obtained exactly. 

In general, the wave function based approach suffers from two major drawbacks: firstly, in some 
situations (as in the previous example) it is impossible to obtain an exact solution of the Schrodinger 
equation which makes it difficult to infer what is the degree of precision of the resulting Wigner 
function. Secondly, even if an exact solution of the Schrodinger equation is known, this method 
only provides the exact Wigner trajectories and in some situations it is desirable that the quantum 
corrections to the classical trajectories might be computed order by order in h. 

The results of this paper provide a new approach to determine the Wigner trajectories for 
bounded systems (satisfying Dirichlet boundary conditions) from an exclusive phase space (or c- 
number) point of view. In fact most of the formalism for the unconfined case can now be applied 
to the bounded case as well. This is certainly true for eq. (69-71,74-76) where one just has to take 
into account that: 1) V(x) is no longer the bulk potential but it also encompasses a boundary 
contribution Vd{x) which is associated with the Dirichlet boundary conditions and was given in 
eq.(53): Vd(x) = ^— {S'(x — b) — S'(x — a)}, 2) in the boundary corrections to eqs. (69,70,74) the 
Wigner function should be evaluated at x + and x~ (in the left and right boundary correction, 
respectively) . 

It is not so clear if the approximative scheme (eq.(77) to eq.(81)) can be also extended to 
the bounded case. This is because the boundary potential deviates strongly from the harmonic 
oscillator potential and yields contributions to the Moyal equation to all orders in Ti. We thus have 
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to worry about the magnitude of these contributions and find a suitable criterion to determine the 
validity of the approximation. These issues are also a problem in the general unconfmed case and 
are certainly not the subject of this paper. Our point here is just that the boundary potential 
approach is, thus far, the unique exact formulation of confined systems that only appeals to phase 
space objects and it thus seems to be the proper framework from which the former or other, more 
suitable, approximative scheme may be devised. 



7.3 Wigner trajectories of a particle in an infinite potential well 

In this section we ignore the issues of validity and convergence of the approximative scheme given by 
eqs. (77-81) and use it to compute the order by order quantum corrections to the classical trajectories 
of a free particle confined to the interval ] — 1, 1[. Quite remarkably, we will see that this (in this 
context, very poor) approximation procedure is already able to reproduce most of the qualitative 
behavior of the exact quantum trajectories. 

The exact Wigner trajectories for a particle in the fundamental state are solutions of eq.(76) with 
the effective potential satisfying (74,75) and were displayed in Fig. 4. The zero order approximation 
to the Moyal equation (53) is given by: 

dF™ P dF™ h 2 fdV R dF™ dV L dF™ \ 

where Vr(x) = S' e (x — 1) and Vl(x) = S' e (x + 1). To simplify the discussion we assume that for 
sufficient small e we have: 

dF^/dp(x-e,p,t) _ dF^/dp(x + e,p,t) _ 
dFF/dp(x,p,t) ~ dFF/dp(x,p,t) ' 

and use these identities to transform eq.(83) into the Liouville equation (77) with V = ^(Yr — Vl)- 
The zero order approximation to the exact quantum trajectories are thus the solutions of eq.(73). 
They are displayed in Fig. 5. where we used the Gaussian approximation to the Dirac delta 
function: S e (x) = e x l e> e' = 0.8e and made e = 0.25 (notice that this approximation was 

used in section 6 to compute the numerical solution of the bounded eigenvalue equation, fig. 2) and 
m — h — 1. For sufficient small e these would be just the classical trajectories. 
We now compute the first order corrections eq.(79). These are the solutions of: 

X= 

m 1 

(84) 

P = £ + !) - S "( x + £ + 1) - ^(x - 1)} , 

where we again ignored the contribution of the factors that involve the (approximate) Wigner 
function. We get the trajectories displayed in fig. 6. 

Finally, we compute the second order corrections (eq.(81)): 

X = 2-, 

P=& + 1) - %(x ~ 1)} + ife {6i 4) (x + 1) - 5^(x - 1)} - ( 85 ) 
l)-5^(x-l)}, 



3840m 



X 
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and get the trajectories displayed in fig. 7. 
Several remarks are now in order: 

a) We see that the first order approximation already displays the deflection of the quantum trajec- 
tories towards regions of higher momentum as the particle approaches the boundaries. 

b) In the second approximation the "islands" appear. They are not confined, as in the exact case, 
to regions of positive or negative momentum, but they are already confined to positive or negative 
positions. Furthermore, the deflection to regions of high momentum is now an exclusive feature (as 
in the exact case) of the higher momentum trajectories. The slower ones are essentially classical. 

c) Not all Wigner trajectories will be confined. For sufficiently high momentum they will escape 
from the box. This is because we are not using the true Dirac delta but a Gaussian approximation 
to it. The free Wigner trajectories correspond to the part of the wave function that is outside from 
the box in fig. 2. 

d) The former approximative scheme can be refined at least in three different directions: firstly by 
taking into account the contributions of the factors that are proportional to the Wigner function. 
Our, rather simple approximative scheme, is in fact, state independent and so could not possibly 
reproduce the exact Wigner trajectories. This step will be more demanding from the numerical 
point of view since one will be asked to solve the partial differential equation for the approximate 
Wigner function before computing each order correction to the classical trajectories. Secondly, by 
increasing the degree of approximation to the Dirac delta function, i.e. by decreasing the spread 
e. This is mandatory if one is interested in the behavior of the quantum trajectories of higher 
momentum. Lastly, and quite obviously by computing higher order corrections in H. 

8 Conclusions 

In this paper we studied the Wigner- Weyl formulation of a particle confined to a finite interval and 
subject to Dirichlet boundary conditions. We found that: 1) the standard procedure of solving 
the unconstrained system and then imposing the boundary conditions does not work out in the 
Wigner- Weyl formulation; 2) the boundary conditions and the confinement of the wave function 
determine a boundary correction to the *-genvalue equation and to the Moyal equation; 3) Dirichlet 
boundary conditions on the confined wavefunction imply both Dirichlet and Neumann boundary 
conditions on the corresponding confined Wigner function. 

Our main task was then the evaluation of the contributions of the boundaries to the *-genvalue 
equation (6) and to the time evolution equation (5) for a wave function satisfying Dirichlet boundary 
conditions. These contributions were shown to have the form of boundary potentials added to 
the Hamiltonian, and the same potentials were obtained in the operator formulation of quantum 
mechanics, when we derived the Schrodinger eigenvalue equation for the confined wave function 
i]){x) = <f>(x)0(x — a)6{b — x). Together with the boundary conditions they are responsible, both 
in the Wigner and in the Schrodinger formulations, for the confinement of the system and the 
consequent discretization of the energy levels. Finally, we extended Baker's converse construction 
to the bounded case, proving the full equivalence between the operator-Schrodinger and the Wigner 
formulation of the confined eigenvalue problem. 

Lastly, we used our previous results to approach the problem of evaluating the Wigner trajec- 
tories of bounded systems and ventured the possibility of applying these methods to non-linear 
systems displaying classical chaotic behavior. 
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The procedure described in this work can be generalized to higher-dimensional multiparticle 
systems and also to any *-genvalue equation of the form A w (x, p)*F^ v (x, p) = aF^(x, p), associated 
with the eigenvalue equation: A(x,p)ip a (x) = aip a (x). From a more general point of view the results 
of this paper provide a first approach to the problem of solving star-equations with boundaries. This 
is a difficult but important problem as such equations are known to play a key part in several fields 
of research ranging from standard topics in quantum mechanics to some important developments 
in M-theory [M. 
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Fig.l: Unconfined fundamental solution of the free particle 
eigenvalue equation satisfying Dirichlet boundary conditions 
at x = ±1 and for h = m — 1 
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Fig. 2: Numerical solution of the first approximation (eq.67) to the 
confined eigenvalue equation: e = 0.25, H — m — 1 and E = n 2 /8 
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Fig.3: Numerical solution of the second approximation (eq.68) to 
confined eigenvalue equation: e = 0.0025, h — m — 1 and E = n 2 
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Fig.5: Zero order or classical approximation to the 
Wigner quantum trajectories: e = 0.25 and H — m — 1. 
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Fig. 6: First order approximation to the Wigner quantum 
trajectories: numerical solution of eq.(84) with e = 0.25 and H — m— 1. 
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Fig. 7: Second order approximation to the Wigner quantum 
trajectories: numerical solution of eq.(85) with e = 0.25 and H — m 
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